Superconducting/ferromagnetic diffusive bilayer with a spin-active interface: a 

numerical study 
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We calculate the density of states (DOS) in a diffusive superconducting/ferromagnetic bilayer with 
a spin-active interface. We use a self-consistent numerical treatment to make a systematic study 
of the effects of the Spin-Dependence of Interfacial Phase Shifts (SDIPS) on the self-consistent 
superconducting gap and the DOS. Strikingly, we find that the SDIPS can induce a double gap 
structure (DGS) in the DOS of the ferromagnet, even when the superconducing layer is much thicker 
than the superconducting coherence lenght. We thus obtain DOS curves which have interesting 
similarities with those of Phys. Rev. Lett. 100, 237002 (2008). 

PACS numbers: 73.23.-b, 74.20.-z, 74.50.+r 
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I. INTRODUCTION 

Superconducting/ferromagnetic (S/F) hybrid struc- 
tures give rise to a fascinating interplay between two an- 
tagonist electronic orders. The ferromagnetic exchange 
field E ex privileges one spin direction while standard su- 
perconductivity favors singlet correlations. This leads to 
a rich behavior which has triggered an intense activity 
in the last years (see e.g. Refs. [l|0). In particular, the 
"superconducting proximity effect", i.e. the propagation 
of the superconducting correlations in ferromagnets, has 
been widely studied. This propagation is accompanied 
by spatial oscillations of the superconducting order pa- 
rameter, because E ex induces an energy shift between 
electrons and holes. As a result, one can build electronic 
devices with new functionalities, such as Josephson junc- 
tions with negative critical currents 3 -, which could find 
applications in the field of superconducting circuits^. 
From a fundamental point of view, it is very instructive 
to study the density of states (DOS) of S/F structures. 
So far, this quantity has been less measured 6 -^"^ than 
critical temperatures or supercurrents. However, this 
way of probing the superconducting proximity effect is 
very interesting because it provides spectroscopic infor- 
mation. One striking consequence of the spatial oscilla- 
tions of the order parameter in the F layer is that the 
zero-energy DOS can become larger than in the normal 
state for certain ferromagnet thicknesses 6 . 

The behavior of S/ F hybrid circuits depends crucially 
on the properties of the interfaces between the S and F 
materials. In this paper, we focus on the case of diffusive 
structures. Diffusive S/F interfaces have been initially 
described with spin-independent boundary conditions 11 . 
It has been found that the amplitude of the supercon- 
ducting proximity effect directly depends on the tunnel 
conductance Gt of an interface (see e.g. Ref. |lj). Later, 
spin-dependent boundary conditions have been intro- 
duced, in the limit of a weakly polarized ferromagne t 12 ' 13 . 
Due to the Spin-Dependence of Interfacial Phase Shifts 
(SDTPS) i 1S ) one has to take into account new conduc- 



tance parameters G^ and G^ at the F and S sides of the 
interface, respectively. It has been shown that G^ and 
G? can significantly affect the behavior of S/F hybrid 
circuits. For instance, G^ can shift the spatial oscilla- 
tions of the superconducting order parameter—. More re- 
cently, it has been found that G^ can induce an effective 

Zeeman splitting A e J^ in a superconducting layer with a 
thickness ds smaller than the superconducting coherence 
lengthscale £s^£. This induces a double gap structure 
(DGS) in the S and F densities of states. However, in 
practice, the regime ds > £s is frequently reached (see 
e.g. Refs. Il9lf20h . Remarkably, DGSs have been recently 
observed at the F side of Ni/Nb bilayers with ds much 
larger than £s>^, although Ref. [H has found that A e J^ 
scales with dg 1 in the low ds regime. Whether a DGS 
persists in the large ds regime is therefore an important 
question, especially in the light of this recent experiment. 

In this paper, we study how G^ and G^ modify the 
DOS of a S/F bilayer. We use a numerical treatment to 
explore a wider parameter range than in previous works. 
In particular, we can reach the limit of thick supercon- 
ductors and larger values of G^. We find that G| shifts 
the spatial oscillations of the superconducting order pa- 
rameter in F, like G^. It can also significantly affect 
the amplitude of the superconducting gap. When ds in- 
creases, the SDIPS-induced DGS becomes narrower, in 
agreement with Ref. 1 181 Nevertheless, it can surprisingly 
persist in the large ds limit. Indeed, on a distance of the 
order of £s near the S/F interface, the resonance energies 
of the S spectrum remain spin-dependent because quan- 
tum interferences make the superconducting correlations 
sensitive to the SDIPS. This behavior is transmitted to 
the whole F layer due to the superconducting proxim- 
ity effect. We thus obtain, at the F side of S/F bilayers, 
DOS curves which have interesting similarities with those 
of Ref. [Tcj, although ds 3> £s . More generally, our results 
could be useful for interpreting experiments. 

This paper is organized as follows. Section II defines 
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the S/F bilayer problem studied in this article. Sec- 
tion III explains the principle of our numerical treat- 
ment. Section IV presents a detailed study of the SDIPS- 
induced DGS. Section V shows the effects of the SDIPS 
on the self-consistent superconducting gap and on the 
oscillations of the zero-energy DOS with the thickness of 
F, Section VI discusses the data of Ref. [l(| Section VII 
concludes. 



II. DESCRIPTION OF THE S/F BILAYER 

We consider a diffusive S/F bilayer consisting of a 
standard BCS superconductor S for —ds < x < 0, and 
a ferromagnet F for < x < (If- We characterize the 
normal quasiparticle excitations and the superconduct- 
ing condensate of pairs with Usadel normal and anoma- 
lous Green's functions Gn,a- = sgn(w n ) cos(# nj(T ) and 
F n .a = sin(#„ !CT ), with 9 n .<j{x) the superconducting pair- 
ing angle, which depends on the spin direction o e {T, J.}, 
the Matsubara frequency w„(T) = (2n+l)7rfcsT, and the 
spatial coordinate a>21. The Usadel equations describing 
the spatial evolution of 9 n ^ write 



and 



i 



2 d 2 n ,a \uJ n \ . . . A(x) 



* s dx 2 ' A 
in S and 



— — sin(6>„ :C 



cos(0„, CT ) (1) 



in F, with 



k n ,a- = V 2 [i<rsgn(uj n ) + (\uj n \ /E ex )\ 



(2) 



(3) 



In the above Eqs., A denotes the bulk gap of the S mate- 
rial, £5 = (fi-F>s/2A B cs) 1 ^ 2 the superconducting coher- 
ence lengthscale, £p = (hDp / 'E ex ) x ^ 2 the magnetic co- 
herence lengthscale, Dp(s) the diffusion constant in F(S) 
and E ex the ferromagnetic exchange field of F. The self- 
consistent superconducting gap A(x) occurring in Eq. ((T|) 
can be expressed as 



A(z)log[^] 



<76{T,i} 

\uj n \<n D 







A(x) 



(4) 



with the Debye frequency of 5, T° = A exp(£)/7rfcs 
the bulk transition temperature of S, ks the Boltz- 
mann constant, T the temperature and £ the Euler con- 
stant. The above equations must be supplemented with 
a description of the boundaries of S and F. We use 
d9 n;Cr /dx\ x= _ d + = d6 n><7 /dx\ x=d - = for the external 

sides of the bilayer. For the S/F interface, we use the 
spin-dependent boundary condition s 12 ' 18 
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dx 



=0 + 



= It sin[0£ CT -0^+i7|Vsgn(w n ) sin[0£j 

(5) 



d9„ 



dx 



x=0- 



11T sin[^ (T -6'^ (T ]-i7^crsgn(w n ) sin[6> 



S 1 



(6) 

with 6 F a = 6 7ha (x = 0+) and 0* ff = 6 n>a (x = 0"). 
These equations involve the reduced conductances 7t = 



Gt£,f I 'Aof and 



G^ S '£f(s)/ a vf(S), the barrier 



asymmetry coefficient 7 = £,s&f/£,fO's, the normal state 
conductivity <Jf(S) °f the F(S) material, and the junc- 
tion area A. Note that we have used a definition of 7^ 
which differs from that of Ref. [HI, to ensure a symmetric 
treatment of 7J and 7^ in Eqs. (0[6|). The microscopic 
expressions of the conductances Gt, G f and G^ can be 
found e.g. in Ref. UH The term Gt corresponds to the 
usual tunnel conductance of the interface. The terms G F 
and G^ can be finite only in case of a Spin-Dependence 
of Interfacial Phase Shifts (SDIPS). The SDIPS results 
from the fact that the scattering phases picked up by 
electrons upon scattering by the S/F interface can de- 
pend on spin due to the ferromagnetic exchange field or 
to a spin-dependent interface potential. Thus, in princi- 
ple, any kind of S/F interface can have a finite SDIPS. 
However, the exact values of G^ and G^ are difficult to 
predict because they depend on the detailed microscopic 
structure of the interface. One possible approach is to 
consider G^ and G^ as fitting parameters which have 
to be determined from proximity effect measurements. 
Note that the derivation of the boundary conditions JH- 
|6|) assumes a weak transmission probability per channel 
(tunnel limit), which seems reasonable considering the 
band structure mismatch between most S and F mate- 
rials. It furthermore assumes that the system is weakly 
polarized. However, there is no fundamental constraint 
on the amplitudes of Gt, G^ and G^ because these pa- 
rameters consist of a sum of contributions from numerous 
conducting channels. 

In S/F circuits, long-range triplet correlations (be- 
tween equal spins) can occur when the circuit in- 
cludes several F electrodes or domains with non-colinear 
magnetizations 22 . Recently, it has been found that this 
effect can also arise in S/F circuits with spin-active inter- 
faces, due to spin-flip interfacial coupling terms which are 
due e.g. to some misaligned local moments at the S/F 
interface 2 ^. In our work, we consider interfaces which 
are "spin-active" in the sense that the SDIPS is finite. 
However, we assume that there is no interfacial spin-flip 
coupling and that F is uniformly polarized. Hence, we 
don't obtain any long-range triplet component with our 
model. 



III. NUMERICAL TREATMENT OF THE 
PROBLEM 

Equations fl][6]) have already been solved numerically 
with a self-consistent procedure in the case 7^ = 7^ = 
(see e.g. Ref. In this paper, we study the case of 
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7^ and 7^ finite, using a numerical treatment based on 
a relaxation method. This treatment is divided into two 
steps. We first calculate the values of A(x) and n ^ 
sclf-consistently with a relaxation method in imaginary 
times. Then, we determine the pairing angle 9 a (e,x) 
corresponding to the calculated A (a;) by using a similar 
relaxation method in real times, i.e. we use w„ = — ie + T 
and sgn(w„) = i in Eqs. (Q][H]), with e the energy, and T = 
0.05A a rate which accounts for inelastic processes 2 ^. 
Finally we obtain the DOS N(s,x) — J^ CT N a (e,x) at 
position x by using N a (e,x) = (No/2) Re[cos[0 CT (s, x)]], 
with Nq/2 the normal DOS per spin direction. Through- 
out this numerical treatment, we use a discretized space, 
with a step of 0.001£;s(f) in S(F). In the following, we 
mainly focus on Np(e) = N(e, x = dp). Ref.[Hhas stud- 
ied analytically S/F bilayers with ds < £s/2, 7^ <C 1, 
7t <C 1, and dp > £f- Our approach allows to go be- 
yond this regime. Note that in Figures Q] to the re- 
sults are shown for E ex — lOOAo, £1_d = QOlksT and 
fc B T = 0.lA . 



IV. THE SDIPS-INDUCED DOUBLE GAP 
STRUCTURE 

A. Variations of the bilayer spectrum with the 
thickness of S 

The left and right top panels of Fig. Q] show the densi- 
ties of states N(e, x — —dg) and N(e, x — 0~) at the left 
and right side of the superconductor respectively, for dif- 
ferent values of ds, and the bottom panel of Fig. Q] shows 
the corresponding DOS Np(e) at the right side of F. For 
ds = 0.35£s, the curves calculated with our numerical 
code (black full lines) are in close agreement with the 
analytical solution given in Ref. [H (black dashed lines). 
The DOS is similar at the two sides of S and displays a 
DOS which reveals the existence of an effective Zeeman 
splitting of the form 

c r s 

l\ z - 2i\ — 7^ - F TH — (7) 

as 

with E TH = HDs I d 2 s the Thouless energy of the S layer 
and Gs — asA/d s its normal state conductance. The 
DGS is also visible in Np(e) due to the proximity ef- 
fect. It becomes narrower when ds increases, in agree- 
ment with Eq. ([7]) which indicates that A^ scales with 
dg . For very large values of ds, the DOS N(e, x = —dg) 
at the left side of S tends to the bulk value Re(cos(#o(e))), 
with 6*o(e) = arctan(Ao/(— is + T)) (see top left panel, 
blue full line). However, a DGS remains clearly visible in 
Np(e), a result which is quite counterintuitive consider- 
ing the low d s expression Eq. (|7|) (see bottom panel, blue 
full line). Note that in the S layer, with the parameters of 
Fig. [TJ ds ^S> £s an d £ = [e = A ], N(s, x) decays from 
its bulk value to N(e, x = 0~) on a distance of the order 
of £s [2£s] near the interface (not shown). In the large 




FIG. 1: Densities of states N(e, x = — d~t> ) (top left panel) and 
N(e,x = 0~) (top right panel) at the left and right side of 
the superconductor respectively, and density of states Nf(e) 
at the right side of the ferromagnet (bottom panel), plotted 
versus e, for different values of ds- The full lines correspond 
to our numericai results. The black dashed fines correspond 
to the anaiyticaf predictions of Ref. [l8l. for ds/^s = 0.35. 



ds limit, the DOS N(e,x — 0~) at the left side of the 
S/F interface does not show a clear DGS for the weak 
value of 7^ used in Fig. [TJ because of the strong DOS 
peak at e = Ao (see top right panel, blue full line) . How- 
ever, a DGS would appear more clearly in N(e, x = CP) 
for larger values of 7?, e.g. using 7^ = 0.4 (not shown). 
The DGS thus seems to persist at large values of ds due 
to an effect which involves a S area with thickness ~ £s 
near the S/F interface and the whole F. 



B. Variations the bilayer spectrum with the 
thickness of F 

Due to the ferromagnetic exchange field E ex , the zero- 
energy DOS Np(s = 0) oscillates around its normal state 
value Nq when dp increase a 6 i 26 ' 27 i 28 . In the large ds 
limit, we find that the DGS can occur at the F side 
for both an ordinary (Np(s = 0) > No) and a reversed 
(Np(e = 0) < A^o) DOS. However, its visibility varies 
with dp, like in the limit ds < Cs/2 of Ref. [18. Fig- 
ure [2 left panel, shows Np versus e, for different values 
of dp and ds/£s = 5. From dp = 0.53£f to 0.58£f, 
Np(e = 0) — A*o is positive and the visibility of the 
DGS increases (see black, blue and green full lines). For 
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e/A Q e /A 

FIG. 2: Density of states Nf(e) at the right side of the ferro- 
magnet, plotted versus e, for different values of cLf (left panel) 
and different values of 7^ (right panel). In the left panel, for 
dp/^F = 1.6, we have multiplied (N F (s) — No) /No by a factor 
50 for visibility of the curve. 



dp = Np(s = 0) — No is negative and the DGS is 
not visible anymore (see red full line). For (If = 1-6£_f, 
the DGS is visible again, with both the inner and outer 
peaks of the DOS inverted due to iV>(e = 0) - iV < 
(sec pink dashed line). 



D. Analytic description of the thick 
superconductor limit 

In order to have a better insight on the persistence of 
the SDIPS-induced DGS at large values of ds, we provide 
in this section an analytic description of the case where S 
is semi-infinite. For simplicity, we assume that the super- 
conducting gap is only weakly affected by the presence 
of the F layer, i.e. A(x) = Ao- We furthermore assume 
that the proximity effect is weak, i.e. B n ^[x e F) « 1 
and 9 n , a (x E S) - 9° < 1, with 0° = arctan(A /|w n |). 
In this case, the Usadel Eqs. fHH]) lead to: 

9 n ,a{x) = 9 na cosh I — I / cosh I ' I 

(8) 

for x E F and 

e n!<7 (x)=0° n +S9l a e W (^ (9) 

for x E S, with r\ n = (l + (uj n / A ) 2 ) ly ' 4 ■ We have intro- 
duced in the above Eqs. 56% a = 9 n ,a{x — 0~) — 0° and 
9^ a = 9 n ,a{x = + ). The linearization of the bound- 
ary conditions (JSJIHI) with respect to these two quantities 
leads to: 



C. Variations of the bilayer spectrum with 

In the parameters range investigated by us (with in 
particular E ex > A, 0.35 < d s /£,s < 10 and 0.4 < 
d F /S,F < 4), no DGS occurs when 7^ = 0. The DGS 
studied in this article thus seems to be a direct conse- 
quence of 7? ^ for ds/^s large as well as ds/£s small. 
Figure [2l right panel, shows the variations of Njr(e) with 
7?, for a constant value of dp. For 7^ = 0, no DGS 
appears. For a very small 7^ (see red full line, corre- 
sponding to 7^ = 0.2 or G^ — Gt), Np(e) shows a 
change of slope which corresponds to a smoothed DGS, 
near e = Ao (from the previous section, this DGS occurs 
only for certain values of dp). When 7^ becomes suffi- 
ciently large, Np{£) shows a clear DGS, i.e. two peaks, 
one above and one below e = Ao, while a local minimum 
is visible for e ~ Ao . The distance between the two peaks 
of Np(e) increases with 7^. In Fig. [21 when 7^ becomes 
too large (7^ > 0.8), the sign of Np(e = 0) — N changes. 
This suggests that 7^ does not only induce DGSs but also 
shifts the oscillations of Np(e) with d F . This last effect 
will be investigated in more details for e — in section 
fVl With the parameters of Fig. [2 the DGS is not visible 
anymore when 7^ becomes larger than approximately 1. 
In the general case, this threshold strongly depends on 
the different parameters characterizing the S/F bilayer. 



and 



9 s 

n.a 



7T(sin(fl T °J + cos(g°)^) 
7t cos(6>°) + i-fF asgn(uj n ) + B n>(7 

77t + Z7|crsgn(w„) 



v; 



77t + i7^crsgn(w„) 



KJ 

An 



(10) 



(11) 



with B nM = fe n>cr taxih[dpk nta /^F]. Importantly, Eqs. |[8}- 
fTT|) are valid provided 59% a <C 1 and 9f la <C 1, which 
requires 7r <C 1, 7^ < 1 and g?f > £ F . We have used 
these hypotheses to simplify Eq. (jTTJ) . The validity of the 
approximation A(x) = Ao will be discussed in section fVl 

Figure[3]shows Np(e) calculated from the analytic con- 
tinuation of Eqs. (fSl TTTj) (black dashed lines) and from our 
numerical code (red full lines), for a weak value of 7, and 
7^ = (left panel) or 7^ =^0 (right panel). The two 
calculations are in relatively good agreement. For the 
parameters used in Fig. [3l we have checked numerically 
that the approximation A(x) — Aq gives results in very 
good agreement with the full resolution of Eqs. ([flS]). At 
e ~ Ao , small discrepancies arise between the predictions 
of the numerical code and of Eqs. (|51 fTT|) . due to resonance 
effects which make 69% a and 9% a larger than for e = or 
e 3> Ao. Equations (|5j|TTj) allow to recover the fact that 
a DGS can appear in Np(e), due to 7^ ^ (right panel). 
In the limit ds S> £s, the pairing angle of the system can- 
not be put under the form 9 a (x, e) = 9(x, e— [oA^/2]), 
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10 12 

e/A e/A 

FIG. 3: Density of states Nf(s) at the right side of the ferro- 
magnet, plotted versus e, for 7^ = (left panel) and 7^ finite 
(right panel). The red lines are calculated from Eqs. (IHl llip 
and the black dotted lines correspond to the self-consistent 
numerical resolution of Eqs. (|H6[I . 

contrarily to what has been found for ds < £s/2. There- 
fore, the notion of SDIPS-induced effective Zeeman split- 
ting is not valid for thick S layers. Nevertheless, from 
Eqs. (|9llip . near the S/F interface, the resonance ener- 
gies of the S spectrum can be spin-dependent because 
quantum interferences make the superconducting corre- 
lations sensitive to the SDIPS on a distance of the order 
of £,s/Vn near the S/F interface. From Eqs. (|8I10|> . this 
behavior is transmitted to the whole F layer due to the 
proximity effect. From Eq. (fTTj) . the energy scale related 
to the occurrence of the DGS has the form: 

A SDIPS = 2A ^ (12) 
G s 

with Gs — crsA/^s the normal state conductance of a 
slab of thickness £s of the 5 material. Interestingly, this 
expression has a form similar to Eq. |J7|). with ds replaced 

by £s (one has 2A = hD s /£? s = E§ H ). Note that for 
a ballistic S/F single channel contact, the SDIPS also 
produces a spin-dependent resonance effect^. However, 
in this case, one does not obtain a DGS but rather a sub- 
gap resonance in the conductance and the zero-frequency 
noise of the system. 

V. SELF-CONSISTENT SUPERCONDUCTING 
GAP AND ZERO-ENERGY DOS OF F 

For completeness, we now discuss the effects of the 
SDIPS on the self-consistent superconducting gap A(x) 
and the zero-energy DOS Np(e — 0) versus dp- 

It is already known that the amplitude of A(x) de- 
creases when "fx or 7 increase, similarly to what happens 
in a S/normal metal bilayer— . Figure 0] compares the 
effects of 7 (left panel) and 7^ (right panel) on A(x) (it 
only shows the effect of 7^ > 0, but the effect of 7^ < 
is similar). One can see that A(x) significantly decreases 
when |7^| increases. Similarly, in a clean superconductor 
connected to a ferromagnetic insulator (FI), A(x) has 



d^ s =5, d^ F =0.8,y T =0.08,^ F =0 




-5 -4 -3 -2 -1 -5 -4 -3 -2 -1 



FIG. 4: Self-consistent superconducting gap A(x) versus the 
spacial coordinate x, for different values of 7 (left panel) and 
J$ (right panel). 

Full[empty] dots: N(e)<[>]N , from numerical calculation 



Full lines: Analytical solution, from Eqs. (8-11) 




FIG. 5: Zero-energy density of states Nf(s = 0) at the right 
side of F, versus the thickness dF of F. The density of states 
calculated numerically is shown with red dots. The full and 
empty dots correspond to Nf(e — 0) < and Nf(e — 0) > 0, 
respectively. The density of states given by Eqs. (|8llll[) is 
shown with black dashed lines. Panel (a) corresponds to a 
case with no SDIPS (7^ = 7^ = 0). Panels (b) and (c) show 
the effect of a finite . Panels (d) and (e) show the effect of 
a finite 7^, in comparison with panel (c) where 7^ = 0. With 
the parameters used here, Eqs. (|8I11|I are in agreement with 
our numerical code only when 7^ = 0. 



been predicted to decrease due to the spin-dependence 
of the reflection phases against FI—. In contrast, in the 
regime of parameters investigated by us, 7? has a negli- 
gible effect on the value of A(x) because it does not occur 
directly in the boundary condition (fTTj) at the S side of 
the interface. From this brief study of A(x), we conclude 
that the approximation A (a;) = Ao used in section HV Dl 
is valid only for sufficiently small values of 7^ , 7 and 7^ . 
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Figure [5] presents the effects of 7^ and 7^ on the vari- 
ations of N F (e = 0) with dp- The DOS calculated nu- 
merically is shown with symbols and the DOS given by 
Eqs. (fSl fTTj) is shown with full lines. In panels (a), (b) 
and (c), we have used 7^ = 0, so that the two calcula- 
tions are in close agreement. In panels (d) and (e), the 
two calculations strongly differ because 7^ is too large 
for the hypotheses leading to Eqs. ([5]fTT|) to be valid. We 
recover the fact that, in the regime d,F > £f, Nf(e = 0) 
shows exponentially damped oscillations with dp&2&. In 
the regime dp < £>, the oscillations of Np(e = 0) with 
d F are less regular. This can be understood from the an- 
alytic description of Section IIVD1 For dp > £p , one has 
B n ,a ~ fcn,cr, so that Nf(s = 0) depends on dp through 
the cosh (df fc„ iCr /^_F) term of Eq. only. For dp < £,f, 
B n a and thus strongly depend on dp, which com- 
plicates the variations of Nf(e = 0) with dF and leads 
to more irregular oscillations. Ref. [l3| has already shown 
that 7^ can shift the oscillations of Np(e = 0) with dp. 
Panels (b) and (c) confirm this result and also shows that 
a positive (negative) 7^ decreases (increases) the ampli- 
tude of Np(e = 0). From panels (d) and (e), 7^ can also 
significantly shift the oscillations of Np(e = 0) with dp, 
in agreement with Fig. [21 right panel. For the param- 
eters used in Fig. [5l 7^ does not modify spectacularly 

the amplitude of Nf(e). For larger values of the 
amplitude of the superconducting proximity effect would 
significantly decrease due to a reduction of A(x) (not 
shown) . 

VI. DISCUSSION ON THE DATA OF PHYS. 
REV. LETT. 100, 237002 (2008) 

We now consider the DOS measurements realized by 
SanGiorgio et al. for Nb/Ni bilayers with d$ = 50 nmi£. 
From Ref. [19| which considers samples fabricated by the 
same team, one has £s ~ 10 nm, so that ds/£s ~ 5. 
However, double gap structures have been observed by 
SanGiorgio et al., which motivates a comparison with 
our model. 

Figure [5] compares the data measured for dp = 1.5 nm 
(black squares) with our numerical calculation (blue 
full line). Our calculation reproduces almost quanti- 
tatively the experimental curve. We have used dp = 
1.5 nm, ds[£s_= 5 and T — 280 mK, in agreement 
with Refs. Il0lil9l . We have also used the exchange field 
E ex = 78 meV, estimated by Ref. [lCj, and the Debye 
temperature T D = 275 K of Nb, taken from Ref. [M We 
have assumed £f = 2.83 nm, -fp = 0.06, 7^ = —1.1, 
7 | = 0.5, and T = 0.025A . Note that from Ref. 0J, 
one has Og 1 ~ 15.9 /xf2.cm and Up 1 ~ 9.7 /iil.cm, so 
that one should have 7 = t;s<Jp /£,f(Ts ~ 2.9 with the 
above values of £s an d £f- In Figure [SI we have used a 
value 7 = 2, which is in relatively good agreement with 
this estimate. The values of 7t, 7J, 7^ and 7 used in 
the fit yield G£/G T ~ 18 and Gf /G T ~ 4. A theoretical 



prediction of these ratios is very difficult because they de- 
pend on the detailed microscopic structure of the Nb/Ni 
interface. However, with a simple delta function barrier 
model, it already is possible to find situations where G^ 
and G^ are larger than Gp (see Appendix A and Fig. 6 
of Ref. [H). Therefore, we think that the SDIPS param- 
eters used by us are possible. 

We cannot reproduce quantitatively the data obtained 
by SanGiorgio et al. for all values of dp with the same 
set of parameters as for dF = 1-5 nm. We think that 
this might be due to the fact that some characteristics 
of the samples like e.g. E ex and thus £f, 7^, 7%, and 
7 can vary with dp^. In the data of SanGiorgio et al., 
from dp = 1.5 nm to dp = 3.0 nm, the distance between 
the two peaks of the DGS increases, like in our model 
(see Fig. [2[ left panel). However, the outer peak of the 
DGS remains very close to e ~ Ao, which seems difficult 
to reproduce with our model. Note that in Ref. [l^, for 
dp = 3.5 nm, the outer peaks of the DOS are inverted, 
whereas a sharp dip occurs at low energies, which can give 
the impression that the innner peaks of the DOS persist 
but are not inverted, in contrast to what we find (see sec- 
tion IV.B). However, we think that the observation of this 
zero-bias dip is not totally reliable. Indeed, SanGiorgio 
et al. explain that, in the DOS of their thickest samples 
(for dp > 3.5 nm), "the zero bias peak is due to the steep 
voltage dependence of the background conductance and 
is therefore a by-product of the data normalization pro- 
cedure" . The zero bias-dip of the dp = 3.5 nm sample 
occurs on the same energy scale as these zero-bias peaks, 
and it goes together with a strange zero-bias singularity 
similar to those shown by the thickest samples. There- 
fore, we are not sure whether a proper interpretation of 
the data of Ref. [13 must take into account this feature. 
Ref. [H suggests that the DGS observed by Sangiorgio 
could be due to triplet correlations. However, it is dif- 
ficult to know whether this interpretation can be more 
satisfying than ours because Ref. [34| does not shown any 
quantitative interpretation of the data and does not dis- 
cuss, for instance, the evolution of the inner and outer 
peak positions with dp. 



VII. CONCLUSION 

In summary, we have calculated the density of states 
(DOS) in a diffusive superconducting/ferromagnetic 
(Sy F) bilayer with a spin-active interface. We have used 
a self-consistent numerical treatment to make a system- 
atic study of the effects of the Spin-Dependence of Inter- 
facial Phase Shifts (SDIPS). We characterize the SDIPS 
with two conductance- like parameters G^ and G^, which 
occur in the boundary conditions describing the S and 
F sides of the interface, respectively. We find that the 
amplitude of A(x) significantly decreases if G^ is too 
strong, whereas it is almost insensitive to G^. In con- 
trast, both Gf and G? can shift the oscillations of the 
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FIG. 6: Comparison between the data of SanGiorgio et al., 
Phys. Rev. Lett. 100, 237002 (2008), for d s = 1.5 nm, and 
our numerical calculation with Ao = 1.3 meV, E ex = 78 meV, 
k B T = 280 mK, d F /£_ F = 0.53, d s /^s = 5, 7 t = 0.06, 7 J = 
-1.1, 7 | = 0.5, 7 = 2, T Debye = 275 K and V = 0.025A . 



zero-energy DOS of F with the thickness of F. Remark- 
ably, we find that the SDIPS can produce a double gap 
structure in the DOS of F, even when the S layer is much 
thicker than the superconducting coherence lenght. This 
leads to DOS curves which have striking similarities with 
those of Ref. More generally, our results could be 
useful for interpreting future experiments on supercon- 
ducting/ferromagnetic diffusive hybrid structures. 
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